On spatial variations of the electron-to-proton mass ratio in the Milky Way 
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We aim at applying the ammonia method, recently proposed to explore the electron-to-proton 
mass ratio ^ = mc/mp, to Galactic sources. Emission lines due to ammonia are originated in the 
numerous interstellar molecular clouds allowing to probe fi at different galactocentric distances. 
High quality radio-astronomical observations of molecular cores in lines of NH3 (J, K) = (1, 1), CCS 
Jn = 2i — lo, HC3N J = 5 — 4, and N2H"'' J = 1 — are used to measure the relative radial 
velocity offsets, AV, between the NH3 (1,1) and other molecular transitions. Robust statistical 
analysis is applied to the n = 207 individual measurements in the Perseus molecular cloud (PC), 
the Pipe Nebula (PN), and the infrared dark clouds (IRDCs) to deduce the mean value of AV and its 
uncertainty. The measured values of AV from the PC, PN, and IRDCs show statistically significant 
positive velocity shifts between the line centers of NH3 and other molecules. The most accurate 
estimates are obtained from carefully selected subsamples of the NH3/CCS pairs observed in the 
PC (n = 21) and the PN (n = 8), AV = 36 ± 7stat ± IS.S^ys m s"^ and 53± ll^tat ± IS.S^ys m s"^ , 
respectively, and from the n = 36 NH3/N2H+ and n = 27 NH3/HC3N pairs observed in the IRDCs, 
AV = 148±32stat ±13.6sys m s~^ and 115±37atat ±31sys m s~^ , respectively. Being interpreted in 
terms of the electron-to-proton mass ratio variation, this gives A/i//i~ (4— 14) x 10~®, which is an 
order of magnitude more sensitive than previous astronomical constraints on this quantity. U Afi/fi 
follows the gradient of the local gravitational potential, as suggested in some scalar field models, then 
the obtained result is in sharp contradiction (~ 5 orders of magnitude) with laboratory atomic clock 
experiments carried out at different points in the gravitational potential of the Sun. However, the 
measured signal might be consistent with chameleon- type scalar field models which predict a strong 
dependence of masses and coupling constants on the ambient matter density. The confirmation of 
a spatial variation of A/i//i will help to distinguish between many theoretical proposals suggested 
to explain the nature of dark energy. New control measurements of the velocity offsets involving 
other molecules and a wider range of objects along with verification of rest frequencies are needed 
to investigate different instrument systematics and to come to a more certain conclusion. 

PACS numbers: 06.20.Jr, 32.30.Bv, 95.36-l-x, 95.85.Bh, 98.38.Bn, 98.38.Dq 



I. INTRODUCTION 



A fundamental characteristic of the universe is that 
decelerating expansion changed to accelerating at red- 
shift 2 ~ 1 [U, [3, Hi ■ This late time acceleration is usu- 
ally attributed to some negative-pressure component of 
the bulk of energy density, referred to as 'dark energy'. 
The nature of dark energy remains unknown in spite of a 
good deal of theoretical and experimental efforts [j| . The 
most popular scenario is a dynamically evolving scalar 
field ip, named 'quintessence' with the energy den- 
sity subdominant at a matter dominated (decelerating) 
epoch (1 < z < 1000), and dominant at lower redshifts 
0. Scalar field models allow to alleviate the so-called 
'cosmic coincidence' problem - the fact that matter and 
vacuum energy contribute by comparable amounts to the 
energy density at the present cosmological epoch Q . The 
scalar field(s) coupling to ordinary matter leads unavoid- 
ably to long-range forces and variability of the physical 
constants violating the Equivalence Principle. Thus, any 
laboratory or space-based experiments as well as astro- 
nomical studies aimed at probing the variations of the 
physical constants, are the most important tools to test 
new theories against observations. 



Present experimental facilities allow us to probe the 
variability of the fine-structure constant a = /{he) and 
the electron-to-proton mass ratio ji = rric/mp, or dif- 
ferent combinations of the proton gyromagnetic ratio Qp 
with a and fj, 

The most accurate laboratory constraint on the tempo- 
ral variation of a of d/a = (— 1.6±2.3) x 10^^^ yr^^ is ob- 
tained from atomic clocks experiments jOj] . Being linearly 
extrapolated to redshift z ~ 2, or i ~ 10"'^'' yr, this labora- 
tory bound leads to Aa/a = (— 1.6±2.3) x 10^^ which is 
well below a conservative upper limit on jAa/al < 10~^ 
stemmed from the analysis of quasar absorption-line sys- 
tems [lS[ni[ll- Here Aa/a = {a'-a)/a, with a, a' de- 
noting the values of the fine-structure constant in the lab- 
oratory and the specific absorption/emission line system 
of a galactic or extragalactic object (the same definition 
is applied to Afj./fi). However, laboratory experiments 
and quasar absorption spectra investigate very different 
time scales and different regions of the universe, and the 
connection between them is somewhat model dependent 

[IlEl. 

The direct and model-free estimate of time variation 
of /i was deduced from comparison of the frequency of a 
rovibrational transition in SFg molecule with the funda- 
mental hyperfine transition in Cs: /i/^ = (3.8 ± 5.6) x 
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10~^^ yr~^ [iBl- Astrophysical estimates of fj, at high 
redshifs are controversial in the claims of nonzero and 
zero changes of /i: A/i//i = (—2.4 ± 0.6) x 10^^ was 
inferred from the analysis of the H2-bearing clouds at 
z = 2.6 and 3.0 but this value was not confirmed 
later [13, llBl- A better accuracy bound was obtained 
at lower redshift z — 0.68 from the analysis of radio- 
frequency transitions in NH3, CO, HCO+, and HCN: 
Afi/fi = (-0.6 ± 1.9) X 10"^ which favors a non- 
variability of n at the level of~3xl0~^. 

In these studies it was assumed that the rate of time 
variations dominates over the rate of possible spatial vari- 
ations. However, this may not be the case if scalar fields 
trace the gravitational field inhomogeneities [13, l2ll |. 
The spatial variation of constants in a changing grav- 
itational potential of the Sun were studied in labora- 
tory measurements of the atomic clock transition ^Sq- 
■^Pq in neutral ^^Sr relative to the Cs standard and of 
the radio-frequency El transition between nearly degen- 
erate, opposite-parity levels of atomic dysprosium (Dy). 
The values for the couplings of a and to the gravita- 
tional field were reported to be ka — (2.3 ± 3.1) x 10^^, 
Lj = (-1.1±1.7) X 10-5 22], and fc„ = (-8.7±6.6) x lO^^^ 
[231 , respectively. 

Thus, the current laboratory and astrophysical esti- 
mates of the temporal and/or spatial changes of the cou- 
pling constants and fermion masses do not exhibit any 
meaningful features which could be interpreted in terms 
of scalar field effects. On the other hand, we do ob- 
serve the late time acceleration of the universe driven 
most probably by ultra-light scalar fields. To recon- 
cile these results chameleon-type models were suggested 
[H m, m, H [28,, M, 30] which allow scalar fields to 
evolve on cosmological time scales today and to have si- 
multaneous strong couplings of order unity to matter, as 
expected from string theory. The key assumption in this 
scenario is that the mass of the scalar field depends on 
the local matter density. This explains why cosmological 
scalar fields, such as quintessence, are not detectable in 
local tests of the Equivalence Principle since we live in 
a dense environment where the mass of the field can be 
sufficiently large and (^-mediated interactions are short- 
ranged. On the other hand, the cosmological matter den- 
sity is about 10"^° times smaller and thus the mass of the 
scalar field can be very low, of order of Hq, the present 
Hubble constant. This allows the field to evolve cosmo- 
logically today. Laboratory experiments aimed at detec- 
tions of chameleon particles are nowadays carried out at 
Fermilab [p. 

In this paper we explore the astrophysical implications 
of the scalar field effect on the electron-to-proton mass 
ratio, fi, following suggestions on experimental tests of 
A/i//i in our Galaxy as formulated in We consider 
the constraints on the spatial variations Afi/fi caused 
by the relative Doppler shifts of molecular transitions in 
cold molecular cores observed in the Milky Way disk at 
different galactocentric distances. 

Our assumptions and basic physical properties of dense 



TABLE L Molecular transitions and Vlsr uncertainties, e„, 
used in calculations of Afi/jj,. The numbers in parentheses 
correspond to la errors (see Sect. II V"^ for more details). 



Transition 




GHz 


Arcst , 

cm 


m s"'^ 


CCS Jjv = 2i — lo 
NH3 {J,K) = (1,1) 
HC3N J = 5 - 4 
N2H+ J = 1 - 


22.344033(1)" 
23.6944955(1)" 
45.4903102(3)' 
93.173777(4)' 


1.34 
1.27 
0.66 
0.32 


13.4 
1.3 
2.8 

13.5 


"Rest frequency from Ref. 
''Rest frequency from Ref. 


M 
H 



molecular cores which will be needed in the sequel are 
discussed in Sect. [Hi The statistical analysis of the line 
position measurements in the Perseus molecular complex, 
the Pipe Nebula, and the infrared dark clouds are pre- 
sented in Sect. IIIIl The obtained results are discussed in 
Sect. IIVI and conclusions are given in Sect. IVl 

II. INTERSTELLAR MOLECULES AS PROBES 
OF Afi/fi SPATIAL VARIATIONS 

Extremely narrow molecular lines observed in cold 
dark clouds provide a sensitive spectroscopic tool to 
study the spatial variations of fundamental constants 
under different environments and gravitational fields. 
Among numerous molecules detected in the interstellar 
medium (ISM), ammonia (NH3) is of particular interest 
due to its high sensitivity to changes in the electron-to- 
proton mass ratio, /x. As was firstly shown for ammonia 
isotopomer ^^NDa [3^ . the inversion frequency in the 
{J,K) — (1,1) level changes as 

^ = 5.6^. (1) 

which is nearly an order of magnitude more sensitive to 
/i- variation than molecular vibrational frequencies which 
scale as £^vib ~ 

The deuterated ammonia shows, however, low abun- 
dances in the ISM molecular clouds, n(ND3)/n(NH3) ~ 
lO^"' — 10~^ [3^. For more abundant ammonia form of 
NH3 , the sensitivity coefficient of the inversion transition 
1/ = 23.69 GHz was calculated in (19i] : 

Au An 

— = 4.46— . 2) 

v ^ 

By comparing the observed inversion frequency of NH3 
(1,1) with a suitable rotational frequency (scaling as 
Etot ~ m) of another molecule arising co-spatially with 
ammonia, a limit on the spatial variation of n can be 
determined p^ : 

^ ^ 0.289 ^'°*"^"^ ^ 0.289— , (3) 
fj, c c 
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AV = V(CCS) - VCNHa) 




Velocity dispersion, ct^CNII^), km s ' 

FIG. 1: Upper panel: CCS (2i - lo) versus NH3(1,1) 
linewidths for cores in the Perseus molecular cloud from the 
total sample of Ref. [33 ]. The error bars represent la stan- 
dard deviations. Lower panel: Subsample of the best quality 
single-component profiles of CCS and NH3 selected from the 
full set of the observed spectra available onhne (see Table . 
In both panels, the dashed and dotted lines are the bound- 
aries for, respectively, pure turbulent, cr„(CCS) — (T„(NIl3), 
and pure thermal, a„(CCS) = 0.55cr„(NH3), line broadening. 
The sample sizes are indicated in the panels. 



where VJ-ot and Vinv are the apparent radial velocities of 
the rotational and inversion transitions, respectively. 

In the present paper we use the comparison of the rel- 
ative radial velocities of ammonia inversion lines and ro- 
tational transitions of another N-bearing (N2H+) and C- 
bearing (CCS, HC3N) molecules (Table HI to set a limit 
on spatial variations of A/i//i. Here we consider some 
physical and kinematic properties of the molecular gas 
clouds where these molecules are detected. 

NH3, N2H+, CCS, and HC3N are usually observed in 
emission in dense molecular clouds (uh > lO** cm~^ ) 
which are represented by a large variety of types in the 
Milky Way disk Two main classes are high-mass 

clumps {Ai > lOOA^o) associated with infrared dark 
clouds (IRDCs), and low-mass cores {A4 < 10A4q) sub- 
divided in turn into protostellar cores containing infrared 
sources, and starless cores (pre-protostellar, or prestellar 
cores) without embedded luminous sources of any mass. 
The latter provide information on the physical condi- 
tions of dense molecular gas preceding gravitational col- 
lapse. Depending on the critical gas density in the core 
center (rina ~ cm~'^), the starless cores can be ei- 
ther dynamically stable against gravitational contraction 
("h2 ^ 10^ cm~'^), or unstable if density exceeds this 
critical value [s^ . The dense cores reside within a larger 
molecular complex or can be isolated. With increasing 
angular resolution the cores reveal a substructure with 
its own range of masses. Being very sensitive to the envi- 
ronmental physical conditions, different molecules probe, 
in general, different substructures within the clouds. 

At high spectral resolution available in the microwave 



FIG. 2: Upper panel: Velocity offset AVccs-nHj versus the 
source radial velocity for points shown in Fig. [TJ lower panel. 
Lower panel: Same as the upper panel but for the points 
marked by asterisks in Table [IT] The vertical error bars in- 
clude both random and systematic errors caused by the un- 
certainties of the adopted rest frequencies (see Table H]). 



range the errors in the molecular line position measure- 
ments are mainly restricted by the uncertainties in labo- 
ratory frequencies, Ei, ^0.1 — 1 kHz, which correspond to 
the Vlsr uncertainties of e,j ~ 1 — 10 m s""'^ [s^l- Taking 
into account that Afi/^^ 0.3AV/c, the expected preci- 
sion of the A/x//x measurements in the ISM is restricted 
by the level of ^ 10~®, i.e. they are about 100 times more 
accurate than the A////Lt estimate deduced at z = 0.68 

This level of accuracy can be achieved, however, only 
ideally since it assumes that molecules are identically (co- 
spatially) distributed within the cloud, and are observed 
simultaneously with the same receiver, beam size, sys- 
tem temperature, and velocity resolution. Violation of 
any of these conditions leads to shifts of the line centers 
which are referred to as the Doppler noise. The input 
of the Doppler noise to a putative A/i//x signal can be 
reduced to some extent if the velocity shifts due to inho- 
mogeneous distribution of molecules and to instrumental 
imperfections are of random nature. 

All molecular transitions listed in Table [J require for 
their excitation high gas density, ~ 10'* cm^'^ , at a 
typical kinetic temperature of Ty^-m ~ 10 K. Ammonia is 
easily detectable in many interstellar clouds. The inver- 
sion transition {J,K) = (1,1) has 18 hyperfine compo- 
nents which can be resolved into nine separate features 
for quiescent nearby dark clouds [40.] . Since hyperfine 
components have different intensities, a single NH3 spec- 
trum allows to determine optical depth, radial velocity, 
intrinsic linewidth, and excitation temperature and thus 
to separate the effects of optical depth and excitation 
temperature 41 1. 

The second N-bearing molecule, which is a good tracer 
of cold and quiescent dense gas, is the molecular ion 
N2H+. The seven hyperfine components of the N2H+ 



4 



TABLE II: Parameters of dense cores in the Perseus molecular cloud [D ~ 260 pc) taken from Ref. [s^l and rep- 
resenting the best quality single component profiles of CCS (2i — lo) and NH3 (1,1). Asterisk marks objects with 
0.55cr„(NH3) < cr^CCS) < cr^NHa) (see Fig.©. 
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J = 1 — transition were studied with high spectral reso- 
lution in [13]. N2H"'" is mainly formed in well-shielded re- 
gions where photoprocesses are unimportant [43| . The in- 
tensity distribution of N2H"'' (1 — 0) in dense cores closely 
matches that of NH3 (1,1), and the relative abundances 
of these molecules are almost constant across the star- 
less core. However, in some cases there were detections 
of ammonia abundances increased by a factor of ~1-10 
towards the core center [i^]- High angular resolution 
observations of nearby prestellar objects show that N- 
bearing molecules trace the inner cores where the density 
approaches 10^ cm~'^ and the carbon-chain molecules 
disappear from the gas-phase because of the freeze-out 
onto dust grains [s^. The observed correlation between 
N2H+ and NH3 velocities and linewidths is another indi- 
cation of probable co-spatial origin of these two species 




The spectral lines of the radical CCS without hyperfine 
splittings and not strongly saturated (r < 3) make CCS 
a well-suited molecule for detailed studies of the velocity 
structure in dark clouds [i^ . CCS lines are rather intense 
in starless, cold and quiescent dark clouds at the earli- 
est stages of molecular cloud evolution, whereas NH3 , as 



mentioned above, tends to be abundant in more chemi- 
cally evolved and dense central regions of molecular cores 
[131. As a rule, CCS shows up as a clumpy distributed 
gas slightly outside these central regions. A pronounced 
spatial anticorrelation in the emission from CCS and NH3 
was observed on large scale rnaps of the Taurus Molec- 
ular Cloud 1 (TMC-1) ridge |48| . pre-protostellar cores 
L1498 [4^ and B68 [50|, as well as on small scale maps 
with angular resolution of ~ 5" {— 1750 AU) in the in- 
frared cloud Bl-IRS where the CCS emission was de- 
tected from three clumps surrounding the central source 
[5l| . The clumps in Bl-IRS exhibit a velocity gradient 
from red- to blueshifted velocities with respect to the 
systemic velocity of the central source which can be in- 
terpreted as interaction of the CCS-emitting gas with 
the molecular outflow. This interaction makes CCS lines 
wider than their thermal widths and comparable with the 
measured widths of the NH3 lines in the Bl-IRS source: 
cr„(CCS) = 0.9 ± 0.1 km s^^ and cr„(NH3) = 1.1 ± 0.1 
km s^"'^ (CCS thermal linewidth is narrower by a factor 
of 0.55 than that of NH3 due to the different masses of 
these molecules). Turbulent broadening could equalize 
the apparent linewidths if both CCS and NH3 trace the 



5 



TABLE III: Parameters of dense cores in the Pipe Nebula 
{D ~ 130 pc) taken from Ref. 65] and representing emission 
line profiles of CCS (2i — lo) and NH3 (1,1). Asterisk marks 
filled circles in Fig. [3] 




FIG. 3: Upper panel: CCS (2i - lo) versus NH3(1,1) 
linewidths for dense cores associated with the Pipe Nebulae 
from the survey of Ref. [6^. No error bars were reported in 
the original paper. The filled circles represent points marked 
by asterisks in Table IIIII The dashed and dotted lines are 
the boundaries for, respectively, pure turbulent, (Ti,(CCS) — 
a„(NH3), and pure thermal, cri,(CCS) = 0.55(t„(NH3), line 
broadening. Lower panel: Velocity offset AVccs-nhs versus 
the source radial velocity set by Vlsh^NHs). 



same volume elements. However, in 5 out of 6 low-mass 
young star-forming regions studied in [s^ the widths of 
the CCS hues are too narrow, cr„(CCS) < 0.55cru(NH3), 
which is consistent with both molecules tracing different 
regions of gas with different kinematics. The only excep- 
tion in this sample is the Bl-IRS cloud where the spatial 
anticorrelation in the CCS- and NHa-emitting gas has 
been revealed [5l|. On the other hand, from the survey 
of dense cores using lines of N2H"'", C3H2, and CCS it was 
deduced [11] that velocities and hnewidths of these three 
molecules are well correlated among themselves and with 
NH3, indicating a kinematic consistency for these species. 

The three strongest hyperfine components of the J = 
5 — 4 rotational transition of HC3N [5J] were partly re- 
solved in radio-astronomical observations [s^. HC3N, as 
well as CCS, is abundant in the early stage of chemi- 
cal evolution of the star-forming regions [56| when the 
ratio of the fractional abundances of HC3N and CCS 
remains almost constant [47| and the spatial distribu- 
tions of these C-bearing molecules match each other quite 
well. The continued chemical evolution leads, however, 
to the spatial variation of the relative molecular abun- 
dances (chemical differentiation) due to adsorption of the 
heavy elements from the gas phase onto grain mantles at 
high densities and dust temperatures < 20 K [sl]. In 
the later stages of the protostellar collapse carbon chain 
molecules are destroyed by high velocity outflows and ra- 
diation from protostars, whereas the same processes favor 
the desorption of ammonia from dust grains [47[ . 

Thus, the spatial distribution of the emitting gas trac- 
ers in star-forming regions is governed by both the phys- 
ical conditions within the dense cores and the time- 
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dependent chemistry. The consistency of the apparent 
linewidths of different molecules is a necessary condition 
for these molecules to be co-spatially distributed, but is 
not a sufficient one. 

Molecular line profiles are also affected by the bulk 
kinematics of the dense clouds. Plane-of-sky bulk mo- 
tions observed in starless cores exhibit velocity gradients 
of ~0.3 to 1.4 km s'^ pc'^ [13, HI]. A differential ro- 
tation and angular momentum evolution in the B68 core 
with velocity gradients of 3.4 and 4.8 km pc~^ for the 
outer and inner material respectively was found in [soj . 
Similarly, the spatial variations of the velocity gradients 
in p restellar and protostellar cores of '^2 km s~^ pc~^ 
[45I ] indicate deviations from simple solid body rotation. 

The gas kinematics along the line-of-sights studied 
in the isolated cores shows that cores move slowly 
( < 0.1 km s~^ ) relatively to their surroundings. Inward 
and/or outward motion with velocity of ~0.1 km s^^ is 
revealed by the skewed line profiles. For instance, blue 
asymmetry of the hyperfine components of the N2II+ (1- 
0) transition seen in two starless cores can be produced 
by absorption along the line-of-sight in the foreground 
inward-moving gas of lower excitation (sij . Acoustic type 
oscillations with amplitude of ~0.1 km s~^ are seen in 
the inward and outward gas motion across the B68 core 
[59| . A redshifted peak of the CCS 82 — 2i emission line 
of ~0.16 km s~^ with respect to the central velocity of 
the NH3 (1,1) line was detected in the LI 551 dark cloud 

The turbulent component of the gas velocity derived in 
[4^ ] from the intrinsic linewidths of the NII3 and N2II+ 
lines is small and subsonic, ctur ^0.1 km s^^ (for Tkin = 
10 K, sound speed is 0.19 km s~^ ). Earlier, a similar 
conclusion has been obtained from NII3 observations of 
low-mass dense cores (6^ . 

The velocity widths in high-mass clumps associated 
with IRDCs are large (~l-3 km s~^ ) and indicate the 
supersonic motion. If the IRDCs are high-mass star- 
forming regions, then, as noted in Ref. [1^, the observed 
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Velocity dispersion, c7^{NHg), krn s ^ 

FIG. 4: Upper panel: N2H+ (1-0) versus NH3(1,1) linewidths 
for massive clumps associated with infrared dark clouds from 
the survey of Ref. [S^ . The error bars represent la standard 
deviations. Lower panel: Same as the upper panel but for 
the points marked by asterisks in Table HVl In both panels, 
the dashed and dotted lines are the boundaries for, respec- 
tively, pure turbulent, cr„(N2lI"'") = o-i,(NIl3), and pure ther- 
mal, 0"„(N2lI+) = 0.77ai, (NII3), line broadening. The sample 
sizes are indicated in the panels. 



large velocity widths are consistent with a model which 
suggests massive star formation in turbulence-supported 
cores with a high accretion rate [1^. The innermost 
regions of such clumps are traced by molecular transi- 
tions with relatively large excitation energies of the upper 
states. To excite such transitions denser and hotter re- 
gions closer to the protostar are required. The NH3 (2,2) 
and (3,3) transitions have, for example, the excitation 
energies of E^'^^ = 65 K and E^-''^ = 125 K, and so their 
velocity distributions should be broader than that of the 
NH3 (1,1) line for which E^'-^^ = 23 K. We may expect 
that other molecules, such as CCS {Jn =1\ — Iq), HC3N 
(J = 5 - 4), and N2H+ ( J = 1 - 0) with excitation en- 
ergies lower than 23 K, would have velocity distributions 
similar to that of NH3 (1,1). 



III. RESULTS 

As noted in the previous section, the observational con- 
straint on A^//L( is set by the velocity offset, AF, mea- 
sured between the ammonia inversion line and a molec- 
ular rotational line, Eq.Q. In this comparison, to min- 
imize possible systematics in the AV^ values, the second 
molecule must follow closely the spatial distribution of 
NH3. The most favorable candidate which fits this con- 
dition is the molecular ion N2H''". As for carbon-chain 
molecules, it seems likely that CCS and HC3N are anti- 
correlated with NH3 , albeit not all of the dense cores dis- 
play chemical differentiation (Ref. Issl. Sect. 4.4). When 
molecules trace the same gas, their turbulent component 
should be equal, and the lighter NH3 should always have 
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FIG. 5: Upper panel: Velocity offset AV'i>j2h+-nh3 versus the 
source radial velocity set by Vlsr(NII3). Lower panel: Same 
as the upper panel but for the points marked by asterisks in 
Table IIVI The vertical error bars include both random and 
systematic errors caused by the uncertainties of the adopted 
rest frequencies (Table H]). 



larger linewidth if turbulent and thermal line broadening 
are comparable, crtur/ctherm ~ 1- If, however, the tur- 
bulent component dominates, the lines may have similar 
widths. We use this criterion to select molecular pairs 
tracing presumably the same gas velocity, i.e. arising 
co-spatially. 

It is to note that the velocity offset AV in Eq.Q is the 
sum of two components, Ay = AV^ + AV^, with AF^ 
being the shift due to /i-variation, and Ay„ the shift in- 
duced by Doppler noise. Assuming the random character 
of the Doppler noise, the signal AV^ can be estimated 
statistically by averaging over a large data sample. The 
validity of this assumption requires, however, additional 
tests as discussed below in Sect. IIII CI 



A. Dense cores in the Perseus molecular cloud 

At first we consider an ammonia spectral atlas of 193 
dense protostellar and prestellar cores of low masses in 
the Perseus molecular cloud [11]. The spectral observa- 
tions of the cores in NH3 (1,1), (2,2), CCS (2i - Iq) and 
CC'^^S (2i — lo) lines were carried out simultaneously 
using the 100-m Robert F. Byrd Green Bank Telescope 
(GET). Each target was observed in a single-pointing, 
frequency-switched mode allowing the entire NH3 (1,1) 
complex to remain within the spectral window. The GET 
beam size at 23 GHz is FWHM = 31" or 0.04 pc at the 
Perseus cloud distance {D ^ 260 pc). The resulting spec- 
tral resolution was 24 m s^^ . Thus, both molecules were 
observed with similar angular and spectral resolutions. 

The typical cores in Perseus have the mean gas density 
~ (1 — 2) X lO'' cm~^ , velocity dispersion cr^ = 0.17 
km s~^ , kinetic temperature Tkin = H K, radius R = 
0.09 pc, and mass M. ^ \ M.q. The Perseus molecular 
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TABLE IV: Parameters of massive clumps associated with infrared dark clouds taken from Ref. 36] and representing emission 
lines of N2H+ (1 - 0) and NH3 (1,1)- Asterisk marks objects with 0.77(t„(NH3) < cr„(N2H+) < cr„(NH3) (see Fig. H lower 
panel) . 





Flsr, km s ^ 






„, km s ^ 


D, 




N2H+ 




NH3 






NH3 


kpc 


G015.05+00.07 MMl 


24.73 


'3 


24.72 


'4; 


2.79(8; 


2.3(1) 


2.5 


G015. 31-00. 16 MM2 


31.12 


*2 


31.05 


3 


0.81(7 


1.34(7 


3.1 


G015. 31-00. 16 MM3 


30.81 


^2 


30.93 


3 


0.62(4 


1.27(7 


3.0 


G019. 27+00. 07 MMl 


26.76 


^1 


26.37 


3 


3.15(3 


2.57(8 


2.3 


G019. 27+00. 07 MM2* 


27.01 


^2 


26.79 


4 


2.62(4 


2.71(9; 


2.3 


G022.35+00.41 MMl* 


52.73 


^2 


52.67 


5 


2.74(5 




3.7 


G022. 35+00.41 MM2 


60.22 


^2 


60.22 


4 


0.83(7 


m 


4.1 


G022. 35+00.41 MM3* 


83.97 


^3 


84.29 


3 


2.23(7 


2.45(7) 5.3 


G023. 60+00. 00 MMl* 


106.89 


^6 


106.10 


5 


3.82(7 


4.1(1) 


6.7 


G023.60+00.00 MM2 


53.34 


^3 


53.49 


3 


3.03(7 


2.47(7) 3.7 


G023. 60+00. 00 MM3* 


105.51 




105.45 


5 


3.4(1) 


3.6(1) 


6.6 


G023.60+00.00 MM4 


53.87 


^5 


53.65 


3 


2.6(1) 


2.20(7 


3.7 


G023.60+00.00 MM7* 


54.00 


^5 


53.74 


3 


1.90(6; 


2.22(7; 


3.7 


G024.08+00.04 MMl* 


113.63 


^1 


113.61 


4 


2.67(3; 


2.9(1) 


7.8 


G024.08+00.04 MM2(1)* 


52.08 


^3 


51.80 


6 


1.5(1) 


1.6(1) 


3.6 


G024.08+00.04 MM2(2)* 


114.45 


^4 


114.02 


5 


2.1(1) 


2.4(1) 


7.8 


G024.08+00.04 MM3* 


51.86 


^6 


51.59 


3 


2.2(2) 


1.89(7) 


3.6 


G024.08+00.04 MM4 


51.95 


^3 


51.59 


3 


0.85(8) 


1.67(7; 


3.6 


G024. 33+00. 11 MMl* 


113.76 


*2 


113.61 


5 


3.53(4 


3.4(1) 


7.7 


G024.33+00.11 MM2* 


118.22 


^4 


118.11 


5 


3.23(9; 


3.5(1) 


7.7 


G024.33+00.11 MM4* 


115.19 


7 


114.89 


7 


3.6(2) 


4.6(2) 


7.7 


G024.33+00.11 MM5* 


117.26 


^4 


117.47 


3 


3.0(8) 


3.06(6) 7.7 


G024.33+00.11 MM6* 


114.43 


^4 


113.85 


5 


2.97(10) 3.2(1) 


7.7 


G024.33+00.11 MM9* 


119.78 


5 


119.28 


3 


1-7(2)^ 


1.90(8) 7.7 


G024. 33+00. 11 MMll* 


112.60 


5 


112.99 


4 


2.18(5 


2.4(1) 


7.7 


G024.60+00.08 MMl 


53.14 


3 


52.99 


3 


3.10(6 


2.37(7) 3.6 


G024.60+00.08 MM2* 


115.16 


^1 


115.11 


4 


2.37(4 


2.3(1) 


7.7 


G025.04-00.20 MMl 


64.00 


^2 


63.53 


2 


2.36(5 


2.08(6 


4.2 


G025.04-00.20 MM2 


63.83 


3 


63.44 


3 


2.77(7 


2.56(8 


4.3 


G025.04-00.20 MM4* 


63.69 


^2 


63.55 


3 


2.18(8 


2.03(6 


4.1 


G034.43+00.24 MMl* 


57.77 


^5 


57.79 


4 


3.15(1 


3.39(9; 


3.5 


G034.43+00.24 MM2 


57.832( 


3) 57.50 


4 


4.04(2 


3.8(1) 


3.5 


G034.43+00.24 MM3 


59.686( 


i) 58.98 


3 


2.90(2 


2.56(7 


3.6 


G034.43+00.24 MM4* 


57.689( 


i) 57.50 


4 


2.81(2 


3.28(9; 


3.5 


G034.43+00.24 MM5 


58.12 


'3 


57.57 


4 


2.58(8 


2.3(1) 


3.5 


G034.43+00.24 MM6* 


58.49 


^2 


58.31 


4 


2.07(7; 


2.59(9 


3.6 


G034.43+00.24 MM8* 


57.63 


5 


57.84 


4 


3.2(1)^ 


3.14(9 


3.5 


G034.43+00.24 MM9* 


59.00 


^2 


58.53 


4 


2.45(6 


2.62(9 


3.6 


118102-1800 MMl 


22.36 


^2 


21.53 


4 


3.49(5 


2.88(8 


2.7 


118151-1208 MMl* 


33.19 


^1 


32.84 


3 


2.10(4 


2.37(8) 


2.8 


118151-1208 MM2* 


29.74 


^2 


29.93 


4 


3.09(4 




2.6 


118151-1208 MM3 


30.65 


^1 


30.86 


6 


0.87(4 


mi 


2.7 


118182-1433 MMl* 


60.00 


^1 


59.72 


4 


3.06(4 


3.18(9 


4.6 


118182-1433 MM2* 


41.10 


3 


40.93 


4 


2.28(8 


2.08(9 


3.6 


118223-1243 MMl 


45.32 


^1 


45.28 


3 


2.16(3 


1.87(7; 


3.6 


118223-1243 MM2* 


45.29 


^2 


45.84 


5 


2.99(6 


2.9(1) 


3.6 


118223-1243 MM3* 


45.73 


^1 


45.66 


4 


2.64(3 


2.95(9 


3.7 


118223-1243 MM4 


45.83 


3 


45.78 


3 


1.97(7 


1.73(8 


3.7 


118306-0835 MMl* 


78.02 


3 


77.88 


3 


2.81(7; 


2.57(8; 


4.9 


118306-0835 MM2* 


76.67 


5 


76.58 


7 




2.3(2) 


4.9 


118306-0835 MM3* 


53.98 


9 


54.45 


6 


ISi, 


2.1(1) 


3.7 


118337-0743 MMl* 


58.37 


^2 


58.37 


5 


3.35(4 


3.7(1) 


3.9 


118227-0743 MM2* 


58.60 


3 


58.53 


7 


2.71(8 


3.2(2) 


3.9 


118227-0743 MM3* 


56.37 


.3 


55.85 




2.44(8; 


3.2(1) 


3.8 



cloud covers an area of 6 by 2 degrees (27 x 9 pc) and 
contains over W^Mq of gas and dust that corresponds 
to the mean gas density of rin ~ 4 x 10^ cm"^ . 

The atlas of 193 sources is very conveniently presented 
on the website (see Ref. HI for details) where the origi- 
nal spectra are shown along with the best fitting models, 
model parameters and their uncertainties. These authors 
found that many cores show evidence for the multiple 
velocity components along the line-of-sight. The total 



number of cores where both CCS (2i — Ig) and NH3 (1,1) 
lines were detected is 98. The central velocities of these 
lines, being averaged with weights inversely proportional 
to the variances of the measurements, reveal that CCS 
lines are systematically offset by 16 m s^^ . It was sug- 
gested in Ref. 3j| that this offset is due to uncertainties 
in the assumed rest frequency of the CCS line. A similar 
marginally significant offset in the velocities of the C^^O 
(1-0) and NH3 (1,1) fines was measured from 18 dark 
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Velocity dispersion, o-„(NH.j), km s" 



FIG. 6: Same as Fig. g] but for HC3N (5 - 4) and NH3(1,1) 
linewidths. Points in tlie lower panel are marked by asterisks 
in Table [V] In both panels the dashed and dotted lines are 
the boundaries for, respectively, pure turbulent, cr„ (HC3N) — 
a„(NH3), and pure thermal, (t„(HC3N) = 0.58o-„(NH3), line 
broadening. The sample sizes are indicated in the panels. 



AV-V(HC^N) V(NH^) 




Source velocity, Vlsr, km s ^ 

FIG. 7: Same as Fig.[5]but for emission lines HC3N and NH3. 
The total sample is shown in Fig. |6l upper panel. The points 
in the lower panel are those marked by asterisks in Table [V] 



cores in which the same position was observed in both 
NH3 (1,1) and C^^O (1-0) [41]. 

Figure m upper panel, shows linewidths of CCS (2i — 
lo) vs. NH3 (1,1 ) for a total sample of data points 
(n = 98) from [3^. In this panel, the dashed and dot- 
ted lines are the boundaries for, respectively, pure turbu- 
lent, cr„ (CCS) = (Tt,(NH3), and pure thermal, (Ti,(CCS) = 
0.55o'.u(NH3), line broadening. A wide scatter of points 
reveals the fact that for many cores a complex struc- 
ture of line profiles hampers the accurate estimate of the 
linewidth. As a next step, we selected from the database 
only those profiles which have been perfectly fitted by 
a single component model. The obtained subsample of 
n = 34 pairs is shown in the lower panel in Fig.[l]with the 
corresponding fitting parameters listed in Table [TTl which 
have been taken from the online database. The origi- 
nal ammonia data were re-reduced with additional small 
Doppler corrections to eliminate IF-to-IF variations (IF 
stands for Intermediate Frequency - a frequency to which 
the radio frequency is shifted as an intermediate step be- 
fore detection in the backend), but there were no signif- 
icant changes from the values derived from the original 
data set [35l| . 

Formal calculations of the weighted means for the to- 
tal sample and subsample give AV ^ — 0.016 ± 0.013 
km s~^ (which is in line with Ref. 1341) and AF™^34 — 
0.040 ± 0.010 km s~^ , respectively. However, to can- 
cel out the Doppler noise component AVn, simple (un- 
weighted) averaging of data is more appropriate (as seen 
from Fig. [21 the AV values have approximately equal 
error bars). The corresponding unweighted values are 
Ay„=98 = 0.044 ±0.013 km s"! and AVn^si = 0.039 ± 
0.010 km s-^ . 

Now let us select only points with fT.u(NH3) > 
cr„(CCS) > 0.55tT.u(NH3), i.e. points for which both emis- 



sion lines arise most probably co-spatially (Fig. (TJ lower 
panel). This gives n = 21 and AV^^^i = 0-050 ± 0.013 
km s-i , and Ay„=2i = 0.048 ± 0.013 km s^^ . 

The distributions of the selected points are shown in 
Fig. [21 Here the vertical error bars include both the error 
of the measurement and the systematic error of = 13.5 
m s~^ , indicating approximate equidispersion of the indi- 
vidual measurements. It is clearly seen that points with 
positive offsets dominate. The scatter of the points re- 
flects complexity of gas kinematics and effects of chemical 
segregation of one molecule with respect to the other in 
dense molecular cores discussed in Sect. [Hi 

The unweighted mean from the n = 21 subsample will 
be considered as the final estimate of the velocity off- 
set for the Perseus dark cores. Accounting for the sys- 
tematic error due to uncertainties in the rest frequencies 
(Table m, we obtain Ay„=2i = 0.048 ± 0.013 ± 0.0135 
km s~^ . The robust redescending M-estimate for the 
shift and scale [H yields AF„=2i = 0.052±0.007±0.0135 
km s^^ (see our final Table fVl)) . 

To understand if there are subclasses among the 
Perseus cores, we also investigated the measured column 
densities of NH3 and CCS. We found the mean ratio of 
iV(NH3)/iV(CCS) = 1180 ±870. There are only two out- 
lying clouds with A^(NH3)/7V(CCS) - 5500 and 25000. 
The plot between the column density ratio versus veloc- 
ity offset does not show any trend. It seems that the dark 
cores are rather homogeneous with respect to the column 
density ratio, and this probably means that there are no 
strong variations in abundances among the clouds. In 
other words there is no evidence for the two different 
populations. 
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TABLE V: Same as Table |IV] but for emission lines of HC3N (5 — 4) and NH3 (1,1). Asterisk marks objects with 
0.58cr„(NH3) < a„(HC3N) < cr„(NH3) (see Fig. El lower panel). 



Source 


Vlsr, km s ^ 




(T„, km s ^ 




D, 




HC3N 


NH3 




rlv^3iN 


NH3 


kpc 


G015.05+00.07 MMl 


24.4(1) 


24.721 


'4' 


] 3.5(3) 


2-3(11 


2.5 


G015.31-00.16 MM2* 


30.92(8' 


1 31.051 


3 


\ l-'llS) 


1.34(7 


1 3.1 


G019. 27+00. 07 MMl 


26.39(5 


1 26.371 


3 


\ 3.7(1) 


2.57(8 


1 2.3 


G019. 27+00. 07 MM2* 


27.02(4* 


1 26.791 


4 


\ 2.9(1) 


2.71(9; 


1 2.3 


G022.35+00.41 MMl 


53.00(7 


1 52.671 


5 


] 3.6(2) 


3.0(1) 


3.7 


G022. 35+00.41 MM3* 


84.36(6 


1 84.291 


3 


\ 1.6(1) 


2.45(7) 5.3 


G023.60+00.00 MMl* 


106.70(6 


1 106. lOl 


5 


1 3.6(1) 


A -X (■\\ 

4.1(1)^ 


D.7 


G023.60+00.00 MM2 


53.37(4 


1 53.491 


3 


\ 3.4(1) 


2.47(7 


1 3.7 


G023.60+00.00 MM4* 


53.91(7 


1 53.651 


3 


\ 2.3(2) 


2.20(7 


\ 3.7 


G023.60+00.00 MM7* 


54.08(5 


1 53.741 


3 


] 2.1(1) 


2.22(7, 


1 3.7 


G024. 08+00. 04 MMl* 


113.68(8 


1 113.611 


4 


\ 2.8(2) 


2.9(n 


7.0 


G024.08+00.04 MM2(1)* 


52.18(7 


1 51.801 


6 


1 0.9(2) 


1 Cil-\\ 

1.6(1) 


3.6 


G024.08+00.04 MM3* 


51.67(4 


1 51.591 


3 


\ 1.6(1) 


1.89(7^ 


1 3.6 


G024.08+00.04 MM4* 


51.90(2 


1 51.591 


3 


\ 1.08(7) 


1.67(7, 


1 3.6 


G024.33+00.11 MMl 


113.69(6; 


1 113.611 


5 


\ 4.0(1) 


3.4(1) 


7.7 


G024.33+00.11 MM2* 


118.1(1) 


118. Ill 


5 


] 3.4(2) 


3-5(l)^ 


7.7 


G024.33+00.11 MM3 


117.2(2) 


116.631 


4 


\ 5-l(4) 


2.16(9 


1 7.7 


G024.33+00.11 MM5 


117.4(1) 


117.471 


3 


] 4.0(3) 


3.06(6 


1 7.7 


G024. 33+00. 11 MM9* 


119.4(1) 


119.281 


3 


\ l-3(2) 


i.yU( 


1 7.7 


G024.60+00.08 MMl 


53.22(5) 52.991 


3 


\ 3.5(1) 


2.37(7, 


1 3.6 


G024.60+00.08 MM2 


115.2(1) 


115.111 


4 


1 3.0(2) 


2.3(1)^ 


7.7 


G025.04-00.20 MMl 


63.67(7; 


1 63.531 


2 


] 3.0(2) 


2.08(6 


1 4.2 


G025.04-00.20 MM2* 


63.63(7 


1 63.441 


3 


\ 2.4(2) 


2.56(8 


1 4 2 


G025.04-00.20 MM4* 


63.76(8 


1 63.551 


3 


\ 2.0(2) 


2.03(6; 


1 4.1 


G034.43+00.24 MMl 


58.04(3 


1 57.791 


4 


\ 3.85(6) 


3.39(9; 


1 3.5 


G034.43+00.24 MM2 


57.86(3 


1 57.501 


4 




3.8(1) 


3.5 


G034.43+00.24 MM3 


60.05(7 


1 58.981 


3 


\ 5.8(2) 


2.56(7; 


1 3.6 


G034.43+00.24 MM4* 


57.50(3 


1 57.501 


4 


\ 3.22(8) 


3.28(9; 


1 3.5 


G034.43+00.24 MM5* 


57.81(8 


1 57.571 


4 


\ 1.8(2) 


2.3(1) 


3.5 


G034.43+00.24 MM6* 


58.34(6; 


1 58.311 


4 


\ 2.5(1) 


2.59(9 


1 3.6 


G034.43+00.24 MM8 


57.7(2) 


57.841 


4 


\ 4.4(4) 


3.14(9 


1 3.5 


G034.43+00.24 MM9* 


58.40(5; 


1 58.531 


4 


\ 2.4(1) 


2.62(9 


1 3.6 


iioiUZ— ioUU IVilVii 


21.53(5 


1 21.531 


4 , 


1 A A ( '\\ 
1 4.4i i 1 


2.88(8 


1 2.7 


118151-1208 MMl* 


33.26(5 


1 32.841 


3 


\ 2.1(1) 


2.37(8; 


1 2.8 


118151-1208 MM2* 


30.30(1 


1 29.931 


4 


] 3.2(3) 


2.9(1) 


2.6 


118182-1433 MMl* 


59.86(3 


1 59.721 


4 


\ 2.58(7) 


3.18(9; 


1 4.6 


118223-1243 MMl* 


45.38(6; 


1 45.281 


3 


] 2.1 1) 


1.87(7, 


1 3.6 


118223-1243 MM2* 


45.5(1) 


45.841 


5 


\ 2.1 3) 


2.9(1) 


3.6 


118223-1243 MM3* 


45.65(7) 45.661 


4 


\ 3.2(2) 


2.95(9; 


1 3.7 


118306-0835 MMl 


78.5(1) 


77.881 


3 


\ 4.0(2) 


2.57(8; 


1 4.9 


118337-0743 MMl* 


58.31(4; 


1 58.371 


5 


1 3.51(8) 


3.7(1) 


3.9 


118337-0743 MM2* 


58.13(6, 


1 58.531 


7 


\ 2.3m 


3.2(2) 


3.9 


118337-0743 MM3* 


55.6(2) 


55.851 


Js; 


\ 3.7(4) 


3.2(1) 


3.8 



B. Dense cores in the Pipe Nebula 

The second sample of molecular lines is based on high 
quality data from Ref. [i^ who observed a starless popu- 
lation of 46 molecular cores in the Pipe Nebula (D ~ 130 
pc). The Pipe Nebula is one of the closest molecular 
clouds with a mass of IQ'^M.q covering an area of 48 
deg^ on the sky. 

The observations were performed with the 100-m 
Robert C. Byrd Green Bank Telescope (GBT). The NH3 
(1,1), (2,2), CCS (2i - lo), and HC5N (9-8) transitions 
were observed simultaneously as in the case of the Perseus 
cloud (the frequency of the HC5N J = 9 — 8 transition 
is 23.96 GHz, and its upper state excitation temperature 
E^^'""^ = 5.8 K). The GBT beam size of - 30" provides the 
angular resolution of ~ 0.02 pc at the Pipe Nebula dis- 
tance, i.e. the internal structure of dense cores (i? ~ 0.09 



pc) was partly resolved. The velocity resolution was 23 
m s~^ . 

There are four sources in the Pipe Nebula sample (Nos. 
12, 40, 101, 109 in Table 4 in Ref. ^ where both NH3 
(1,1) and HC5N (9,8) hnes are observed. The molecule 
HC5N is 4.4 times heavier than NH3, i.e. its lines should 
be intrinsically narrower than ammonia lines if molecules 
trace the same gas. However, each of these sources ex- 
hibits HC5N velocity dispersion broader than NH3: in 
average cr„(HC5N) = (1.5±0.1) (Ti,(NH3). An even worse 
relationship is seen between HC5N and CCS velocity dis- 
persions: in average cr^(HC5N) = (1.8 ± 0.5) cr„(CCS). 
This could be expected if NH3 and CCS trace not all 
the gas along the line of sight but only a part associated 
with the dense compact core. However, not all of the 12 
NH3 and CCS emission sources listed in Table Hill show 
the expected correlation of the velocity dispersion. Four 
sources marked by open circles in Fig. [31 upper panel. 
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FIG. 8: Upper panel: HC3N (5 - 4) versus N2H+ (1 - 0) 
linewidths for massive clumps associated with infrared dark 
clouds [3^ . The error bars represent la standard deviations. 
Lower panel: Same as the upper panel but for the points 
distributed between the dashed and dotted lines which con- 
strain respectively pure turbulent, (t„(HC3N) = a„(N2H"''), 
and pure thermal, (j„(HC3N) = 0.58(t„(N2H+), line broaden- 
ing. The sample sizes are indicated in the panels. 



are those where two species trace different regions. The 
distribution of the remaining 8 points is shown in Fig. [31 
lower panel. Again, as for the sources from the Perseus 
molecular cloud, we find a positive offset in the radial 
velocities. The difference between the CCS and NH3 
velocity has the mean of AVn=s ~ 0.044 km s~^ and 
the standard deviation of 0.020 km s~^ . Adding the 
error of the rest frequencies, we obtain the final values 
AVn=8 = 0.044 ± 0.020 ± 0.0135 km s-^ , and the robust 
M-estimate AVn=s = 0.069 ± 0.011 ± 0.0135 km s"^ 
(Table IVj). 



C. Massive clumps associated with Infrared Dark 
Clouds 

The third example deals with more distant objects — 
cold, dense, and massive clumps {M > IOOMq) of 1-10 
pc sizes associated with the IRDCs which are located at 
distances larger than 2 kpc and seen in silhouette against 
the bright diffuse mid-infrared emission of the Galactic 
plane. The gas kinetic temperature of IRDCs is as low as 
the temperature of the dense molecular cores, Tkin < 20 
K [b^. This means that massive clumps in IRDCs are 
either starless or in a very early stage of high-mass star 
formation. 

A recent survey of 55 massive clumps in NH3 (1,1), 
(2,2), (3,3), N2H+ (1-0), and HC3N (5-4) lines has been 
carried out at the Nobeyama Radio Observatory 45-m 
telescope The N2H+ and HC3N (5-4) lines (Tabled 
were observed simultaneously using two SIS receivers. 
The third receiver HEMT was used for NH3 lines. The 
data were obtained within one day for each position, and 
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Source velocity, km ss ' 

FIG. 9: Upper panel: Velocity offset AyjjQ^j^_f^2jj+ versus the 
source radial velocity set by Vlsr(N2H+). Lower panel: Same 
as the upper panel but for the points from the lower panel 
m Fig. [8] The vertical error bars include both random and 
systematic errors caused by the uncertainties of the adopted 
rest frequencies (Table HI . 

different receivers used in these observations allow the 
Doppler tracking to be completely independent between 
these molecular lines [13 • 

The angular and velocity resolutions in these obser- 
vations were different. At the frequencies of the NH3, 
HC3N, and N2H"'" transitions the half-power beam width 
was 73", 37", and 18", and the velocity resolution was 
0.50, 0.25, and 0.12 km s~^ , respectively. At a distance of 
~5 kpc the spatial structures of ~1.8, 0.9, and 0.4 pc can 
be resolved in, correspondingly, NH3, HC3N, and N2H+ 
lines. The inhomogeneity of the observational parameters 
and significant turbulence in IRDC cores (atmb ^1 — 3 
km s^^ ) increase the dispersion of the Doppler noise in 
massive clumps as compared with low-mass molecular 
cores. 

Table HV] fists 54 sources from ^ where both N2H+ 
(1-0) and NH3 (1,1) fines were detected. These molecules 
are expected to have similar spatial distributions. How- 
ever, Fig. m upper panel, shows that some of the N2H+ 
lines have either too large or too small velocity dis- 
persion as compared to NH3, which may be related to 
a large angular resolution difference (18" vs. 73") in 
observations of these molecules. Such points were ex- 
cluded from the sample, the remaining n = 36 points 
are presented in Fig. 01 lower panel. The velocity off- 
set distributions from the total sample and the subsam- 
ple are shown in Fig. [5l upper and lower panels, respec- 
tively. The weighted and unweighted means are as fol- 
lows: Ay™^3g = 0.124 ± 0.043 km s'^ , and AF„=36 = 
0.122 ± 0.049 km s^^ . In this case the laboratory er- 
ror is equal to = 13.6 m s^^ . For the final esti- 
mate we take the more appropriate unweighted mean 
AVn=36 = 0.122 ±0.049 ±0.0136 km s^^ , and the robust 
Af-estimate AVn=36 = 0.160 ± 0.032 ± 0.0136 km s^^ 
(Table Ell- 
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TABLE VI: Sample mean values (unweighted) AV, standard deviations (Jmis, robust A/-estimates of the sample mean AVm, 
and the scale am deduced from the original data and from the subsamples of molecular lines showing self-consistent linewidths. 
Col. 9 lists velocity offsets corrected for uncertainties of the rest frequencies (see Sect. IIV"Xl for more details). The errors of the 
mean arms/^/n and omI \f'ti are listed, respectively, in columns (4) and (6). 





Molecular 


Sample 


ATT" 








S^sys, A V 


A T /"corr 

AVm , 


Object 


pair 


size, n 


km s 


1 —1 

km s 


km s 


1 — 1 
km s 


1 —1 

km s 


1 —1 

km s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Perseus 


NH3/CCS 


98 


0.044 ± 0.013 


0.129 


0.040 ± 0.007 


0.073 


0.0135 


0.024 




34 


0.039 ±0.010 


0.058 


0.045 ± 0.007 


0.043 


0.0135 


0.029 






21 


0.048 ± 0.013 


0.060 


0.052 ± 0.007 


0.032 


0.0135 


0.036 


Pipe 


NH3/CCS 


12 


0.087 ± 0.039 


0.135 


0.039 ± 0.023 


0.081 


0.0135 


0.023 


8 


0.044 ± 0.020 


0.057 


0.069 ±0.011 


0.030 


0.0135 


0.053 


IRDCs 


NH3/N2H+ 


54 


0.157 ±0.040 


0.294 


0.160 ±0.030 


0.220 


0.0136 


0.148 




36 


0.122 ±0.049 


0.294 


0.160 ±0.032 


0.190 


0.0136 


0.148 




NH3/HC3N 


43 


0.138 ±0.043 


0.282 


0.110 ±0.032 


0.210 


0.0031 


0.105 




27 


0.105 ±0.045 


0.234 


0.120 ±0.037 


0.190 


0.0031 


0.115 




N2H+/HC3N 


41 


-0.056 ± 0.047 


0.301 


-0.020 ± 0.037 


0.240 


0.0138 


-0.013 




22 


-0.033 ± 0.055 


0.258 


-0.017 ±0.034 


0.160 


0.0138 


-0.010 



The correlation between velocity dispersions of the 
HC3N (5-4) and NH3 (1,1) emission lines is shown in 
Fig. [HI In this case the results were obtained with more 
comparable angular resolutions (37" vs. 73"). The total 
sample of the HC3N/NH3 pairs (n = 43) is given in Ta- 
ble El Several of the measured HC3N linewidths demon- 
strate significant supersonic turbulence of the emitting 
gas with cri,(HC3N) > o'i,(NH3). Such lines were omit- 
ted from the following analysis, the remaining n — 27 
points are shown in Fig. [51 lower panel (the selected 
points are marked by asterisks in Table fVl first column). 
Fig. [71 illustrates the distribution of the measured off- 
sets. The calculated statistics are as follows: AV^^27 = 
0.147 ± 0.039 km s'^ , and Ay„=27 = 0.105 ± 0.045 
km . Again, taking the unweighted estimate and 
accounting for the systematic error of £„ = 3.1 m s^^ 
we come to the final mean value of the velocity offset 
Ay„=27 = 0.105±0.045±0.0031 km s^^ , and the robust 
M-estimate AF„=27 = 0.120 ± 0.037 ± 0.0031 km s'^ 
(Table ED. 

The data from Ref. [36] allow us to answer the ques- 
tion whether the revealed positive velocity offsets are re- 
lated to ammonia lines or not. Following the same pro- 
cedure as above, we can compare the velocities obtained 
from the sources which were observed in both N2H+ and 
HC3N. The total sample of n = 41 pairs is shown in 
Fig. m upper panel. The heavier molecule HC3N should 
have linewidths smaller or equal to that of N2H+. How- 
ever, a considerable number of points ('--'50%) lie above 
the dashed line constraining pure turbulent line broad- 
ening. This is expected since N2H+ molecule traces 
more compact central regions of the dense clouds. The 
n = 22 points which presumably trace the same gas 
are shown in Fig. [51 lower panel. The corresponding 
velocity offset distributions are shown in Fig. |9l The 
mean values are AV^^22 — —0.053 ± 0.044 km s~^ , and 
^Vn=22 = — 0.033 ±0.055 km s~^ , and thus the final esti- 
mates are AV„=22 = -0.033 ±0.055±0.0138km s"! , and 



the robust M-value AF„=22 = -0.017 ± 0.034 ± 0.0138 
kms-i (Table Ei- 

In this case both weighted and unweighted values are 
consistent with zero velocity offset. Since N2H+ and NII3 
share similar physical conditions in their formation and 
destruction processes and, hence, trace each other in spa- 
tial distributions, the obtained result indicates that am- 
monia is a key element responsible for the revealed ve- 
locity offset. 

To demonstrate that the positive velocity offset does 
not result from our selection criterion we summarize all 
measurements complemented by robust redescending M- 
estimates for the shift and scale in Table ED Both results 
are in good agreement within the la uncertainty interval. 
The sample mean values reproduce a marginally signif- 
icant signal ( > 3(t) independently on the sample size. 
When ammonia lines are excluded from the analysis, the 
signal disappears. 



IV. DISCUSSION 

In the previous section we showed that different molec- 
ular clouds demonstrate statistically significant velocity 
offset between the ammonia inversion transition (J, K) = 
(1, 1) and rotational transitions of other molecules. Here 
we consider possible reasons of this offset. 



A. Uncertainties of the rest frequencies 

In cold dark clouds with the kinetic temperature Tkin ^ 
lOK the thermal line broadening is very small (uth ~ 50 
m s~"^ for, e.g., CCS). This means that the detailed study 
of velocity offsets requires accuracy of rest frequencies of 
about 1 m s~^ [B^I- Among molecules used in the present 
work only NH3 fulfills this requirement. The hyperfine 
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TABLE VII: Hyperfine components of the HC3N J = 5 - 4 transition. For MoUer (Ref. [O), Yamamoto (Ref. HI), and JPL 
frequencies only the numbers after the decimal point are shown, the integer part is as in Lapinov's data [631 ■ Col. 10 and 11 
list velocity offsets between Yamamoto and Lapinov (Y-L) and between JPL and Lapinov (JPL-L) values. 





Transition 










Rest Frequency (MHz) 




At; 


(m s-i) 


J' 


F' 


J 


F 


Weight 


Lapinov 




MoUer 


Yamamoto 


JPL 


Y-L 


JPL-L 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


(9) 


(10) 


(11) 


5 


5 


4 


5 


0.034 


45488.83643 


'42; 


.8368 


;io) 










5 


4 


4 


3 


0.658 


45490.26138 


42 


.2614 


5) 


.263(1) 


.2580(3) 


10.7 


-22.3 


5 


5 


4 


4 


0.812 


45490.31378 


42 


.3137 


5) 


.316(1) 


.3102(3) 


14.6 


-23.6 


5 


6 


4 


5 


1.000 


45490.33740 


42 


.3373 


6) 


.341(1) 


.3336(4) 


23.7 


-25.0 


5 


4 


4 


5 


0.0003 


45490.63123 


43 


.6317 


16) 










5 


4 


4 


4 


0.034 


45492.10858 


.42; 


.1085 


'9) 











structure of the observed (J, K) ~ (1,1) transition of 
NH3 is known with the accuracy of ~ 0.6 m s~^ [6§|. The 
available laboratory uncertainties for other molecules are 
significantly larger. Besides, there are shifts between rest 
frequencies reported in different publications. Below we 
estimate such systematics and their effect on the mean 
values of the velocity offsets obtained in Sect. IIIII 

The frequency 22.344033 GHz of the CCS Jn = 2i-1o 
transition (Table [l| used in Refs. [3l,[6H| was calculated 
from a comparison with the three strongest hyperfine 
components of the HC3N J = 5 — 4 transition observed 
towards a cold dark cloud L1498 assuming that the spa- 
tial distribution of CCS is the same as that of HC3N 
[5^ . According to Ref. the T^lsr value obtained 

from the HC3N transition was used to determine the fre- 
quencies of the CCS transitions without any corrections, 
and the frequencies of the CCS transitions were calcu- 
lated from Gaussian fitting of the observed profiles. The 
error of 1 kHz indicated in Table H] corresponds to one 
standard deviation obtained in the fit. However, up to 
now the J = 5 — 4 transition has not been measured 
in laboratory. The analysis of all available laboratory 
data on HC3N carried out in Refs. [zO, [zH and inde- 
pendently in Ref. [ssj shows systematic frequency shifts 
(Table lVlI| with respect to the data presented in Ref. [55| 
and in the JPL Catalog (http://spec.jpl.nasa/gov/). On 
the other hand both Ref. 68] and Ref. 71] frequencies 
are in good agreement with the former having slightly 
better accuracy (e^ ~ 2.8 m s~^ ). The mean value for 
the shift between rest frequencies used in Ref. [13] and 
calculated in Ref. fesj is AT4hift = 16 m s'^ (Col. 10, 
Table Em ■ Thus, the presented in Table IVll robust M- 
estimates of the velocity offsets between CCS and NH3 
towards the Perseus molecular cloud and the Pipe Neb- 
ula corrected for this systematic shift will be AT^°2i = 
At/„=2i - AV^hift = 0.036 ± 0.007 ± 0.0135 km s'^ , and 
AV^^^Z'i = AK=8 - AT/shift = 0.053 ± 0.011 ± 0.0135 
km s-i , respectively (Table Ell Col. 9). 

Now let us consider the IRDCs from Ref. [s^. In this 
measurements the lower spectral resolution does not al- 
low to resolve hyperfine structure of the HC3N J = 5 — 4 
transition. The spectra of the HC3N lines were fitted to a 
single Gaussian function and the corresponding radial ve- 



locities were calculated on base of a single hyperfine line 
J = 5-4,F = 5- 4[33] with the frequency 45490.3102 
MHz taken from the JPL Catalog (TableEIIl Col. 9). As 
a result, a small systematic shift due to unresolved hy- 
perfine components was introduced. Its value can be esti- 
mated if we consider the center of gravity for the weighted 
central hyperfine components of the J = 5 — 4 transi- 
tion listed in Table IVlTl The frequency of the center of 
gravity is z^c.g. = 45490.3058 MHz and thus AVc.g. = 29 
m s~^ . However, as seen from Table EUl (Col. 9), the 
JPL fre quen cies are systematically shifted with respect 
to Ref. [6^] data and the mean value of the shift is 
AKhift = — 23.6 m s~^ (Col. 11). Combining both shifts, 
we obtain a correction of 5.4 m s~^ . This means that the 
velocity offset between HC3N and NH3 listed in Table EJ 
should be as follows: AV^^^'^j = 0.115 ± 0.037 ± 0.0031 



km s-i (Table ED Col. 9) 



The rest frequencies for the second molecule N2H+ 
observed towards the IRDCs Ref. [s^ were taken from 
Ref. [12] ■ The hyperfine frequencies of the N2H+ J — 
1 — transition were estimated from observations of 
a molecular cloud L1512 in Taurus in both N2H+ (93 
GHz) and C3H2 (85 GHz) lines assuming that the C3H2 
and the N2H"'" have exactly the same radial velocity. 
The absolute frequencies of the N2H"'" hyperfine tran- 
sitions measured in Ref. [4^ are known with precision 
of 7 kHz. However, a factor of two more accurate 
data are given in the Cologne Database for Molecular 
Spectroscopy (CDMS) described in Re fs. 0,[ll|- Both 
datasets are compared in Table IVIIIl The error of 4 
kHz {sv ~ 13.5 m s^^ ) reported in Table [J corresponds 
to the CDMS. The mean shift between the data from 
Ref. mi and the CDMS frequencies is Al^hift = 12 
m s^^ . Thus, the velocity offset between N2H+ and 
NH3 from Table EI should be corrected to AV^l'L = 



AVr 



n— 36 " 



n— 36 

-0.012 = 0.148±0.032±0.0136 kms^^ , and be- 
tween HC3N and N2H+ to AV^g'i2 = AK=22 + 0.00 66 = 
-0.010 ±0.034 ±0.0138 km s-i (TableEH Col. 9). 



13 



TABLE VIII: Hyperfine components of the N2H''' J = 1 — transition. For Caselli (Ref. 42) frequencies only the numbers after 
the decimal point are shown, the integer part is as in the MoUer data (Ref. iTlh . Col. 10 lists velocity offsets between Caselli 
and MoUer (C-M) values. 







Transition 








Rest Frequency (MHz)t 


A?; (m s~^ ) 


J' 


Fi 


F' J 


Fi 


F 


Weight 


MoUer Caselli 


C-M 


(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 



0.143 
0.604 
0.111 
0.218 
0.159 
0.051 
0.604 
0.111 
1.000 
0.278 
0.134 
0.017 
0.252 
0.100 
0.076 



93171. 
93171. 
93171. 
93172. 
93172. 
93172. 
93173. 
93173. 
93173. 
93173. 
93173. 
93173. 
93176. 
93176. 
93176. 



6171(41) 
9129(42) 
9129(42) 
0499(42) 
0499(42) 
0499(42) 
4748(42) 
4748(42) 
7722(41) 
9657(43) 
9657(43) 
9657(43) 
2611(43) 
2611(43) 
2611(43) 



.6210(70 
.9168(70 

.0533(70) 



.4796(70) 

.7767(70 
.9666(70 



.2650(70) 



12.5 
12.5 

10.9 



15.4 

14.5 
2.9 



12.5 



B. Systematics due to kinematic shifts 

Spectroscopic measurement of the absolute radial ve- 
locities at the level of 50 — 100 m s~^ is subject to many 
physical effects [t^I- Fortunately, relative radial veloci- 
ties are less sensitive to most of them. However, some 
small systematic shifts, in addition to those discussed in 
Sect, ini may be present in the measured values of ISV . 

First of all, the measurement of the relative radial ve- 
locities between two molecules requires the use of a com- 
mon reference frame. In radio-astronomical observations 
it is usually the Local Standard of Rest (LSR). The trans- 
formation of the apparent Doppler velocity to Vlsr in- 
troduces an error of ~ 1 — 3 m s~^ , stemming from the 
following corrections: the observer's motion due to (1) 
the Sun's motion with respect to the LSR, (2) revolution 
of the Earth-Moon barycenter round the Sun, (3) motion 
of the Earth center round the Earth-Moon barycenter, 
and (4) rotation of the Earth /73J . This error ought to be 
added to the total error budget along with the systematic 
errors listed in Table IVIl ^th column if molecular lines 
were not observed simultaneously. 

The second error arises from gravitational shifts of the 
spectral lines due to the annual variation of the observer's 
distance from the Sun because of the eccentricity (e ~ 
0.01671) of the Earth's orbit (A = 1 AU). This produces 
the peak-to-peak variation of the radial velocity (in units 
of the speed of hght, c) of 



5V^ = 



2eGMQ 



(4) 



The numerical value of 5Vg is about 10 cm s ^ , and in 
our case this error can be neglected. 

If we further define the dimcnsionless gravitational po- 
tential of the Sun on the Earth's orbit as (negative sign 



is omitted) 

$0 ^ GMq/{Ac^) + GM(sl{R®c^) ^ 1-06 X 10"* , (5) 

then the gravitational redshift for photons escaping an 
object of mass M and radius R scales as 



/ M 



$0 



(6) 



Different gravitational redshifts between molecular 
lines may occur from the chemical stratification in dense 
molecular cores. However, this effect is negligible because 
of a relatively low gas density and rather large linear sizes 
of the clouds. For a typical starless core shown in Fig. 3 
in Ref. [3^ , the gravitational velocity shifts of NH3 and 
CCS lines are ~ 0.2 cm s^^ and ~ 0.1 cm s^^ , respec- 
tively. 

Systematically redshifted CCS 2i — Iq lines with re- 
spect to the NH3 (1,1) central velocity may occur in con- 
tracting clouds when the near side of the external shell 
of the cloud with the CCS emission is moving towards 
the higher density core as suggested for, e.g, the L1551 
dark cloud [5l|. However, as discussed in Sect. [Hi in 
dark clouds both kinds of kinematic shifts (redward and 
blueward) are observed. Such bulk motions should ran- 
domize the Doppler shifts and provide a zero offset for 
a large size sample of dark clouds. Besides, all molec- 
ular cores in the Perseus molecular cloud are probably 
optically thin in the CCS 2i — Iq emission since CC'^'^S 
was not detected along any of the lines-of-sight [s^l and 
32 s : s = 22 : 1. This means that the CCS emission 
line profile is formed in the whole volume of the cloud 
including both its near and distant sides. This is not the 
case, however, for the NH3 emission. The total opacity in 
the (1,1) line transition varies between 0.84 and 10.9 for 
the cores listed in Table [Til The optically thick clouds 
may cause some systematics of both signs (positive or 
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FIG. 10: Schematic location of the Perseus molecular cloud 
{I ~ 158°, h ~ -20°, D ~ 260 pc), the Pipe Nebula {I ~ 
356°, 6 ~ 7°, D ~ 130 pc), and the IRDCs (£ ~ 15° - 30°, 
& < 1°, _D > 2 kpc) in projection onto the Galactic plane. The 
Galactic center longitude is ^ = 0°, by definition. Note that 
the distance to the IRDCs is not scaled and they lay much 
further out. 

negative depending on the character of the internal bulk 
motions within the core) in the central positions of the 
NH3 lines. However, three cores with r{i.i) < 1 (Nos. 20, 
112, and 151) demonstrate positive velocity offsets with 
the mean value of ^100 m s"-'^ . 

Among analyzed clouds the most accurate measure- 
ments of Ay belong to the Perseus cores. We may use 
the M-estimate from the n = 21 sample (Sect. IIV Ap 
as a tentative value for the found velocity offset, AV ~ 
36 ± 7stat ± 13.5sy, m s"! , or A^^ = 36 ± 15 m s^^ {la 
C.L.). Thus, the considered systematics do not explain 
the revealed offset. Assuming that this offset is caused 
entirely by the change of /i, one obtains from Eq.Q 
A^/^ = (3.5±1.4)xl0-^or |A^/^| < 5 x 10"^ This is 
the most accurate estimate of A/x/// based on the spec- 
tral analysis of astronomical objects. Next we consider 
whether this value is in line with other constraints. 



C. Gravitational-potential dependence of jj, 

First we explore the dependence of the electron and 
proton masses on the strength of the gravitational po- 
tential. Theoretical models dealing with variations of 
physical constants do not require the cosmological tem- 
poral evolution of a and/or to be dominated over their 
local spatial changes in the Milky Way. For instance, 
theories like that of Ref. consider direct coupling be- 
tween the local matter density and the scalar fields that 
drive varying constants. The spatial change of the con- 
stants in these models is proportional to the local grav- 
itational potential and the rate of spatial variation may 
exceed the rate of time variation (see, e.g., Ref. HO and 
references cited therein) . The effect of scalar fields on /i 



is additive (similar to gravity) and can be expressed as 
(Refs.llllTi): 

— = A:^A$ , (7) 

where is a dimensionless coupling constant of a mass- 
less (very light) scalar field to the local gravitational po- 
tential <&, and A$ is the difference of the gravitational 
potentials between two measurement points. 

This class of models suggests that the coupling of the 
scalar field ip to baryonic matter and fields is not too 
strong since otherwise it would produce a scalar-mediated 
'fifth force' that is not observed [2^. It is also assumed 
that the (/s-mediated force is a long-range one. 

The respective contribution from the Galactic gravita- 
tional potential $g {r, C) and the molecular cloud gravita- 
tional potential $cioud to the local gravitational potential 
is the sum 

<i>(r,C) = $G(r,C) + *cioud , (8) 

where r and C denote two of the cylindrical coordinates. 

The gravitational potential of IRDCs (A4 ~ lOOA^©, 
i? 1 pc) and starless cores {M ~ IQMq, Q.l pc) is 
much smaller than that of the Sun on the Earth's orbit, 
$0, namely, <I>cioud ^ 1/2000 $o- The value of ^oi'^X) 
can be calculated from analytic models of the Milky Way 
[73,[2l,[2l| which describe fairly well the galactic rotation 
curves. The galactic coordinates of the analyzed clouds 
are shown in Fig. (TUl The distance of the Sun from the 
center of the Milky Way is ro = 7.94 ±0.42 kpc [80|. The 
Perseus molecular cloud is located towards the galactic 
anti-center, whereas the Pipe Nebula and the IRDCs in 
the direction of the galactic center. Between 3 and 8 kpc 
the galactic potential is about 2 x 10~^, and $0 changes 
within this interval very slowly: A$g ^ 10~^. 

With Afi/fj.— 3.5 X 10~^, one obtains fc^ ^ 1, which 
is in sharp contradiction with estimates based on atomic 
clock experiments: — (-1.1 ± 1.7) x 10~^ [22|. There 
are two possible explanations for such a large discrep- 
ancy: 

1. The measured velocity offset is mainly due to kine- 
matic shifts. Then, taking the atomic clock limit 

< 3 X 10~^ and the change of the galactic grav- 
itational potential A$ 10~^, one obtains from 
Eq. ([7]) an upper limit on |A/z//i| < 10"^^. Unfor- 
tunately, the level of 10~^^ is below the precision 
of modern astronomical observations of molecular 
lines. 

2. The electron-to-proton mass ratio does not linearly 
correlate with metric variations, i.e., k^ = in 
Eq. ([7]). The nonzero velocity offset is then caused 
by the change of the electron and proton masses. 
The corresponding theoretical background for this 
item is discussed below. 
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D. Dependence of ^ on matter density 

A class of models where baryon masses and coupling 
constants are strongly dependent on the local matter den- 
sity, p, and where spatial variations of ^{p) and a(p) 
can be more pronounced than the temporal cosmological 
variations was considered in a number of publications dis- 
cussed in Ref. [30|]- In this approach it is suggested that 
the coupling of a scalar field to matter is much stronger 
than gravitational, but the fifth force problem does not 
arise since the linear coupling of the scalar field to mat- 
ter is effectively suppressed by the matter density itself. 
The dynamics of the scalar field in these models depends 
on p: in low-density environments it is determined by 
the scalar field potential V{ip)^ whereas in high-density 
regions (like the Earth's surface) it is set by the matter-iy9 
coupling. The range of the scalar-mediated force for the 
terrestrial matter densities is then very short, less than 
1 mm 

Because of the strong (^-coupling to matter the mea- 
surements of the frequency drifts in the atomic clock ex- 
periments ought to be insensitive to the changes in the 
gravitational potential at Earth caused by the eccentric- 
ity of Earth's orbit. The minimum of V{lp) would show 
no shifts and Aa/a (Ap/p) would be zero with good 
accuracy. However, comparison of atomic clocks at high 
altitude satellite orbits with laboratory clocks may yield 
useful constraints on model parameters as discussed in 
[30l | who also suggested to test spatial variations of the 
physical constants through atomic and molecular line ob- 
servations of the interstellar clouds in our galaxy. In this 
case the difference of matter densities between the terres- 
trial environment (p® ~ 3 x 10^** GeV cm~'^ ) and, e.g., 
dense molecular cores (jOcioud 3 x 10*^ GeV cm^"^ ) may 
cause a shift in the expectation value of the scalar field 
iprn leading to a change in masses and coupling constants. 

If we interpret the nonzero offset of AV ~ 40 m s~^ 
in terms of the change of fermion masses, then for low 
density environments, such as the interstellar clouds, the 
change in the electron-to-proton mass ratio is given by 

^-(ee-ep)f. (9) 

In this model the scalar field is coupled to the standard 
model fermions through the dimensionless constants 
and ^p, with being normalized to one. 

The factor (^c — 1) can be much larger than unity in 
accord with astrophysical constraints discussed in 30] . If 
we assume that (^c — l)/2 > 1, then </j^j < 3.5 x 10^*, or 
'Pm<2x 10-4. 

Thus, considered observations of the interstellar molec- 
ular clouds being interpreted in terms of the Ap/p vari- 
ation are in agreement with models which treat the ip- 
mediated force as a short-range force depending on the 
matter density. In accord with model predictions, the 
estimated value of the scalar field is much less than 



one. The measured spatial variation of Ap/p does not 
contradict the laboratory studies on atomic clocks (see 
Sect. lIV C|) since extremely different density environment 
in the terrestrial measurements and in the interstellar 
medium prohibits direct comparison of the obtained re- 
sults. 

There are several extragalactic measurements of Ap/p 
in which gas densities similar to those in the Milky Way 
clouds are observed. These results, being interpreted in 
terms of the spatial variations of Ap/p , can be compared 
with our estimate of Ap/p = (3.5 ±1.4) x 10"^: it is 
consistent with the limit on \Ap/p\ < 3 x 10"^ derived 
from the NHa-bearing clouds at z — 0.68 but in 
conflict with the value Ap/p = (— 24±6) x 10~^ obtained 
from molecular hydrogen II2 absorption lines at z 3 
p^ . However, the same H2-bearing clouds at 2 3 
as analyzed in Ref. jT^ were studied in Ref. [17]. The 
variability of p was not confirmed at the level of \ Ap/ p\ < 
4.9 X 10"^. Recently a more stringent limit on Ap/p was 
found at z ~ 3 in Ref. [Q: Ap/p = (2.6 ± 3.0) x 10"^. 
Since at densities as low as observed in the interstellar 
medium the dependence of p on p is extremely weak, 
the value of Ap/p in quasar absorbers is expected to 
be at the same level as in the interstellar clouds, i.e. ~ 
10^^ — 10^ This level of accuracy can be probed with 
the next generation of high resolution spectrographs such 
as ESPRESSO planned for the Very Large Telescope U. 

E. Temporal vs spatial variability of p 

Our measurements of Ap/p being formally divided by 
the distance to the sources in light years imply a time 
variation of Ap/p of about 10^^^ — lO^^'' yr^^. This 
variation rate is orders of magnitude larger than the lab- 
oratory bound quoted in Sect. HI \p/p\ < 10"^'^ yr~^, and 
than the extragalactic bounds based on the upper limit 
on Ap/p at z ~ 3, \p/p\ < 6 x 10"^*^ yr~^ [3, and 
at 2 = 0.68, \p/p\ < 4 x 10-^^ yr'^ [l|. Such glar- 
ing discrepancy may serve as another argument in favor 
of chameleon-type models: due to principal dependence 
of couplings and masses on matter density the labora- 
tory bounds and astronomical measurements cannot be 
directly compared. Taken at face values, we may sug- 
gest that the rate of the spatial variations of p exceeds 
significantly the rate of temporal variations. 



V. CONCLUSIONS 

In this paper we propose to use the ammonia inversion 
transitions in conjunction with low-lying rotational tran- 
sitions of other molecules to probe the spatial changes of 
the electron-to-proton mass ratio p from observations of 
molecular clouds located in the disk of the Milky Way at 
different galactocentric distances. 

The reported estimates of Ap/p are obtained for the 
Perseus molecular cores |34i, the Pipe Nebula [1^, and 



16 



the infrared dark clouds [s^ . We found evidence for the 
spatial variability of fj, at the level of A/x//i ~ (4 — 14) x 
10~® which is the most accurate astrophysical value to 
date. The statistical reliability of this result seems to be 
high enough but some instrument systematics could be 
possible. The revealed positive velocity offset between 
the NH3 inversion transition and rotational transitions 
of other molecules should be confirmed by independent 
observations at different radio telescopes. 

The result obtained can be compared with other tests 
of spatial variations of physical constants in the solar sys- 
tem which are based on atomic clock laboratory measure- 
ments. Thus, if A/i//i follows the gradient of the local 
gravitational potential, as suggested in some scalar field 
models, then our estimate of A/i//i contradicts signifi- 
cantly (~ 5 orders of magnitude) the atomic clock con- 
straints obtained at different points in the gravitational 
potential of the Sun on the Earth's orbit. However, the 
measured signal of A/i//i ~ (4—14) x 10~^ is in agreement 
with chameleon-type scalar field models which predict a 
strong dependence of masses and coupling constants on 
the ambient matter density. 

To be completely confident that the derived veloc- 
ity shift is not due to kinematic effects in the clouds 
but is the appearance of the spatial variation of A/^//i , 
new high precision radio-astronomical observations are 
needed for a wider range of objects. Such observations 
should include essentially optically thin lines of molecules 
co-spatially distributed with ammonia in order to reduce 
the Doppler noise. It is also desirable to increase the ac- 



curacy of the rest frequencies of the CCS (2i — Iq) and 
N2H+ (1 - 0) transitions since their present uncertain- 
ties translate into the radial velocity error of '--^ 13 m s^^ 
whereas the error of the velocity shift between, for in- 
stance, CCS and NH3 in the Perseus cloud can be as 
smah as 5 m s^i (Table [Til)- 

In addition, the search for spatial variations of the fine- 
structure constant a in the Milky Way using mid- and 
far-infrared fine-structure transitions in atoms and ions 
[S^l, or the search for spatial variations of the combina- 
tion of a^//i using the [Cii] A158 /im line and CO rota- 
tional lines [HI would be of great importance for cross- 
checking the results. 
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